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We analyze the t — t'—t" — J model at low doping 8 <C 1 by chiral perturbation theory and show 
that the (1,0) spiral state is stabilized by the presence of t',t" above critical values around 0.2 J, 
assuming t/J = 3.1. We find that the (magnon mediated) hole-hole interactions have an important 
effect on the region of charge stability in the space of parameters t', t" , generally increasing stability, 
while the stability in the magnetic sector is guaranteed by the presence of spin quantum fluctuations 
(order from disorder effect). These conclusions are based on perturbative analysis performed up to 
two loops, with very good convergence. 



I. INTRODUCTION 
AND SUMMARY OF BASIC NOTATION 



The nature of the ground state of doped quantum an- 
tiferromagnets is a central issue for the theory of corre- 
lated electrons^ While it is well established that arbi- 
trary small doping destroys conventional Neel order, it 
is far less clear what is the structure of the emerging 
ground state, and how it co-exists with superconductiv- 
ity. One of the early proposals made in Ref. 0, and later 
explored in the context of the Hubbard and the t — J 
models) 3 ! 4 ! 5 ! 6 ! 7 ! 8 ! 9 ! 10 ! 11 ! 12 ! 13 was that for small doping the 
collinear Neel order gives way to a non-collinear spiral 
state. Energy is gained since the holes can hop easier 
in a spiral background. However it was also established 
that the spiral state is itself unstable with respect to long 
wavelength density fluctuations, i.e. has a tendency to- 
wards phase separation, signaled by a negative charge 
compressibility. The general problem of what states are 
susceptible to instabilities, such as phase separation, and 
what is the true ground state of the doped antiferromag- 
net, is enjoying a lot of current attentionii2*i 4 *i 5 rfi£iiLi2ii 9 . 
Some works favor, under certain conditions, inhomoge- 
neous ground states, such as stripes, where the holes seg- 
regate into ordered structures ^iii whereas others argue 
that ultimately the ground state is uniform. 17-18 

We have recently revisited the problem of stability of 
the spiral state in the extended t—t'—t!' — J model ji 3 - and 
have found that the uniform spiral state is stable (at low 
doping) above certain critical values of t', t" . In addition 
we have shown that superconductivity is supported in the 
stable region. Partial motivation for this research was 
provided by the considerable body of evidence that (in- 
commensurate) magnetism co-exists with superconduc- 
tivity in the LSCO family*^ Such experiments are usu- 
ally interpreted in the context of stripesjSi however the 
(uniform) spiral scenario is quite consistent with the ma- 
jority of data, especially in the range of doping where the 
charge order is "fluctuating" , i.e. not in the ground state. 



Disorder is also expected to be important in these com- 
pounds, and can be readily incorporated into the spiral 
stated 

In the present work we use chiral perturbation 
theorjiSiiSi as our main technical tool. The theory pro- 
vides a rigorous perturbative treatment of long-wave- 
length dynamics with Goldstone quasiparticles in a sys- 
tem with strong interactions. The starting point is the 
ground state of the Heisenberg model [t — J model at 
half filling) which incorporates all spin quantum fluctu- 
ations. The chiral perturbation theory allows a regular 
calculation of all physical quantities in the leading order 
approximation in powers of doping 8. Subleading powers 
of 5 depend on the short-range dynamics and hence can- 
not be calculated without uncontrolled approximations. 
Therefore we cannot determine reliably what is the exact 
value of 5 so that it is small enough to justify our calcu- 
lations. However the limit i5 < 1 justifies the approach 
paramctrically. 

In our previous worki£ we concluded that the spiral 
state is stable in a certain region of parameters t' and 
t" . It was shown that the magnetic stability is due to 
the existence of spin quantum fluctuations (order from 
disorder effect). While the strong renormalization of one 
hole properties due to scattering with magnons was taken 
into account in Ref. 

El 

the calculation of the compress- 
ibility at finite doping (charge response) was made for 
free holes (zero- loop approximation), and the calculation 
of the magnetic response was performed in the one-loop 
approximation. In the present paper we extend our cal- 
culations up to two loops for both the charge and the 
magnetic response. The results demonstrate that pertur- 
bation theory converges very well. The main conclusion 
that the magnetic stability is provided by spin quantum 
fluctuations (order from disorder effect) remains valid in 
two loops, however the higher corrections are found to 
have an important effect on the charge stability region in 
the space of parameters t' , t" , generally increasing stabil- 
ity. 

We start by summarizing some results and the notation 
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of Ref. 13 which will be needed for our calculations. In 
the S=l/2 extended t — J model, one allows hopping 
to nearest- neighbor sites (t), as well as (diagonal) next 
nearest- neighbors tf, and next next nearest sites t" on a 
2D square lattice: 



" — ' ^] C icr C JO 1 C i<J C h° ' ' ,n-' J-''' 
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(1) 



The cJ CT operators act in the space with no douple elec- 
tron occupancy. The number of added holes, per site, 
is denoted by 8, and is referred to as doping. We also 
set the total number of sites N = 1 and measure all the 
energies from now on in units of J which we set to one, 
while t is set to its physical value: 



t = 3.1 



J= 1. 



(2) 



The same-sublattice hoppings t',t" will be considered as 
parameters, having in mind that in the cuprates their 
values are argued to be t' ps —0.8, t" w 0.7, from fits of 
ARPES measurementsi 25 ' 26 

The one-hole Green's function is calculated in the self- 
consistent Born approximation (SCBA), which produces 
reliable results due to the absence of low-order vertex 
corrections i!2*2L24 The hole dispersion minima are lo- 
cated at the points ko = (±7r/2, ±7r/2), which are the 
centers of the four faces of the magnetic Brillouin zone 
(MBZ). In the vicinity of these points the dispersion is 
quadratic: 



Ih 
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(3) 



where k is defined with respect to ko, and k\ is perpen- 
dicular to the face of the MBZ, while &2 is parallel to it. 
The quasiparticle residue near ko is ss Zk a = Z. By 
implementing the SCBA equations numerically, we have 
obtained the coefficients in the following expansion, valid 
in the range -1< t' < 0, < t" < lM 

/3i = 1.96 + 1.15i'+0.06i' 2 + 2.70i" + 0.53i" 2 + 0.50i'i", 
lh = 0.30-1.33i'-0.19i' 2 + 2.80i" + 1.06i" 2 -0.14i'i", 
Z = 0.29 + 0.055i' + 0.195*". (4) 

Then the renormalized hole-magnon interaction is writ- 
ten in the long- wavelength limit, by introducing the op- 
erator 7r q = a q — 0l„i where a q ,/3 q are the two spin 
waves in a two-sublattice antiferromagnet. Calling the 
two sublattices "a" and "b" and introducing the corre- 
sponding (renormalized) hole operators /ika,^k6, we ob- 
tain the Hamiltonian describing the interactions between 
holes and spin waves: 



h.sw 



= Vm, 



k,q 



( /l k+ q( A&7r q + h.c.) 



(5) 



For simplicity from now on we refer to the renormalized 
(dressed) holes simply as holes. The vertex in (JSJ is 
is perpendicular to the MBZ face): 
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(6) 



We stress that M q — ► at q — * in accordance with 
the general Adler's relation^ In a spiral state the spins 
deviate from the Neel order, leading to lowering of the 
ground state energy and creation of a gap in the hole 
spectrum. The spiral pitch Q can be directed along the 
(1,1) or (1,0) direction of the lattice and is proportional 
to doping Q ~ 8. The new hole operators %p, tp, and the 
corresponding energies are: 
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(7) 



The spiral is a coplanar state with spins lying in a fixed 
plane. The arbitrary phase y, is related to the orientation 
of this plane, which, as we will see below, does not appear 
in any physical observables. The operators ifa and <pk 
describe spinless fermions, and due to (jZJ the f fermion 
band is completely empty. The gap and the Fermi energy 
are: 



2ZH 



2+2 



1 



= t^ttV^ 8. (8) 



A, 



Here p s = Z p /A, Z p = 0.72 (S=l/2) is the renormal- 
ized spin stiffness of the undoped antiferromagnet, and 
N p stands for the number of full pockets enclosed by the 
Fermi surface. For the (1,0) spiral N p = 2, while for the 
45° spiral (1,1) there is only one pocket N p — 1,14 No- 
tice that except for the band splitting JjJ which absorbs 
the soft mode of the doped Heisenberg antiferromagnet 
we assume that the dispersion Q remains rigid under 
doping. 

The rest of the paper is organized as follows. In Sec- 
tion II we analyze the charge stability of the spiral states, 
while Section III is devoted to the stability in the mag- 
netic sector. The dispersion variation due to doping is 
discussed in Section IV. Section V contains our conclu- 



II. CHARGE RESPONSE AND CHARGE 
STABILITY 

We now proceed to calculate the charge compressibil- 
ity x, defined as = d 2 E/d8 2 , where E is the ground 
state energy. The two simplest contributions to the en- 
ergy come from the Fermi motion of holes and from the 
energy gain due to the spiral formation. These two con- 
tributions are»14 
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FIG. 1: (a) First order interaction energy. Dashed line rep- 
resents the vertex HUH , (b) Second order direct energy, (c.) 
Second order exchange. The labels ip and ip show which type 
of fermion is present in the intermediate state. 



(a) 




(b) 



FIG. 2: Third order ladder diagrams, (a) Direct, (b.) Ex- 
change. 



Introducing the rescaled variables: p\ — p\\J (5\j {2ep), 
P2 = P2 y//3 2 /(2eF), we have explicitly: 



where the normalized distribution is: 
and we define the function: 



(12) 

;i-iki), 



The superscript (0) in the energy shows that interaction 
effects have not been taken into account (no loops). We 
note that the terms that are not quadratic in S are not ex- 
plicitly written, i.e. we have subtracted from the energy 
the constant and linear terms: Eaf + const. 5. 

The holes interact via spin wave exchange. Treating 
the hole-magnon interaction ||5jl in the second order of 
perturbation theory (the h a , hb operators have to be ex- 
pressed via the ip, ip operators from Q) we obtain the 
effective low-energy hole-hole interaction vertex: 



M 2 

UJ a 



-8Z 2 t 2 ^. 

|q| 2 



(10) 



Here we have neglected both the retardation effects as 
well as the hole energies in the denominators. This ap- 
proximation should work well in the limit of low doping 
S <C 1 since typical hole energies ~ ep ~ 5, while the 
magnon energy w q ~ |q| ~ kp ~ VS ^> S, for 6 <C 1. 
Notice also that Zt w 1 and consequently the interaction 
is strong even in the long- wave- length limit. 

The one-loop correction to the energy is given by the 
exchange diagram in Fig. 1(a). The corresponding ex- 
pression reads 



(11) 



p,k 



where the summations are within one hole pocket en- 
closed by the elliptic Fermi surface, = 6{ep — £k)- 
The direct term, proportional to r q= o is set to zero, by 
keeping in mind the presence of a small non-zero fre- 
quency in the denominator of H10|) . Notice the overall 
sign (+) in Eq. 1|11|) which is opposite to the sign appear- 
ing for an electron gas with Coulomb interactions. This 
is due to the fact that the interaction (|10fl is attractive. 



(/?i//? 2 )fc 2 2 



(13) 



An explicit analytical evaluation of is not possible, 
but the numerical integration in l|12|) does not cause any 
problems. 

The two-loop corrections to the energy are given by 
the diagrams in Fig. l(b,c), which contain both ip — ip 
and ip — tp scatterings. The first type of contribution 
produces: 



E 



(2) 

tp — 1p 



2 ^ ft 

P,k,q P 



r k 



c q c p+k c q+k 

xn p n q (l - %>+k)(l - n-q+k), 
while the second one is: 

i ^ r2-r w r r 



E 



(2) _ 



k^ p+k+q 



2 ^ ft 
P,k,q € P 
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6 p+k ~~ 6 q+k 



(14) 



(15) 



Here both the direct and exchange terms are added up, 
as shown in Fig. l(b,c). In the rescaled variables we have 
the form convenient for numerical integration: 



E 



(2) 

tjj—lp 
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x(l 



d 2 kd 2 pd 2 q 



E 



(2) 
ip-tp 



~ "p+k 
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)(l-^+k)7r7r 



[k.(k + p + q) 



d 2 kd 2 pd 2 q 



q+k 

— /k/ k+p+q 



l [k.(k + p + q) + A/e f ] 



(16) 



Notice that the energy corrections Ijl6(l do not contain 
infrared or ultraviolet divergences in contrast to the 3D 
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FIG. 3: The inverse compressibility x j as denned 
in Eq. lfT%jl . calculated for t' < 0. 



electron gas with Coulomb interaction^ In fact each of 
the two-loop diagrams in Fig. ^ contains an ultraviolet 
logarithmic divergence, In A, where A ~ -4= 3> 1 is an 
ultraviolet cut-off. However the divergences are canceled 
out between the direct and the exchange diagrams, since 
the fermions ip, tp are spinless. 

In the next, third order of perturbation theory the only 
diagrams containing ultraviolet divergences are the "lad- 
der" ones, shown in Fig. [21 Both diagrams separately are 
proportional to In 2 A, and, once again, the divergences 
are canceled out exactly when the direct and exchange 
terms are added together. By examining also higher or- 
ders we find quite generally that the perturbative expan- 
sion behaves well. The absence of divergences is due to 
the dimensionality of the problem as well as the Gold- 
stone nature of the mediator of the interaction. 

We have verified that the hole-magnon vertex only acts 
within a single pocket, i.e. holes from different pockets 
do not interact. Thus, the total ground state energy in- 
cluding one and two-loop corrections is: 



E = 25<o> + N^+Ef^+Ef^) 



(17) 



All terms in this expression are proportional to S 2 and all 
the terms scale as 1/N p since expressions (|12|l and 
contain e 2 F . Therefore, defining x 1 to be the inverse 



compressibility of the (1,0) spiral state, 
ix _1 5 2 , we have: 



En 



(1,0) 



E, 



(14) 



2E, 



(i,o) - X 



,-1x2 



(18) 



Thus whenever the (1,0) state is lower in energy -E(i,o) < 
-E(i.i), it is stable x _1 > 0, whereas the (1,1) state is 
always unstable, because if < ^(i.o) thenE 1 ^ ^ < 

and hence Xn i) < 0- 

In Eqs. H17I18I) we have neglected the contact inter- 
action energy E c , arising from the last (density-density) 
term in proportional to J = 1. In fact one could 
imagine adding an additional nearest neighbor Coulomb 
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FIG. 4: Region of stability of the spiral state in the t—t'—t" — 
J model at t = 3.1 (throughout the work we take J = 1). The 
unstable region, below the critical line, is marked as x" 1 < 0. 
Inset: Stability boundaries calculated in different orders of 
perturbation theory. The dot-dashed line shows the result 
obtained in Ref. [ljj without the interaction corrections (zero 
loops). The dashed line is the stability boundary calculated 
in the one- loop approximation (E^ and E^ in Eq. 1171 ). 
The solid line, which is the same in the inset and in the main 
figure, shows the stability line calculated in the two-loop ap- 
proximation. 



repulsion energy V > to of the form Hy = 
VJ2(ij) n i n j- Then the total contact energy to lowest 



order is: Ef 11} = 0, Ef 1Q) = (V - l/2)Z 2 S 2 . The pres- 



ence of the V term generally increases x _1 ; however as 
long as V ~ 1, the contact energy's contribution to the 
total energy is extremely small, since all the terms in 
Eq. I|17|) contain powers of the dominant scale t ps 3. 

We have found numerically that the relative impor- 
tance of each successive order of perturbation theory de- 
creases roughly by a factor of 3. For example, we estimate 
the first order with respect to the kinetic (Fermi motion) 
energy: N P E^/E F w 1/3. The second-order contri- 
bution relative to the first order is also approximately: 



\E 



(2) 

tp—tp 



E^'_ v \/E^ fa 1/3. We have therefore stopped 



at the second order of perturbation theory and expect 
the contribution of the next orders to be negligible. 

The inverse compressibility is shown in Fig. 3 as 
a function of t' for two values of t" . When the point 
X -1 = is reached, the system becomes unstable. The 
stability analysis leads us to the phase diagram of Fig. 4: 
the spiral phase is stable above the solid line. In Fig. 4 (in- 
set) we show the boundaries of stability calculated in dif- 
ferent orders of perturbation theory. The dot-dashed line 
is the zero loop result, taking into account E^ ' only (no 
interaction effects) , while the dashed line is the boundary 
calculated in the one-loop approximation (E^ and E^ 
in Eq. (|17pV Finally, the solid line is the stability bor- 
der calculated in the two-loop approximation. The insert 
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demonstrates the good convergence of the perturbation 
theory. It is interesting that account of the interaction 
corrections (one and two-loop corrections) enlarges the 
region of stability in the space of parameters t',t". 



III. MAGNON STABILITY 

We now proceed to verify that the region which is sta- 
ble with respect to density fluctuations (x _1 > 0), shown 
in Fig. 4, is also stable in the magnetic sector. The start- 
ing point is the spin-wave Green's function*^ 



D n (u,q) = - 



2ljJq [u 



2 Wq P„( W ,q)] 2 -4c^|P a ( w ,q)| 2 ' 



(19) 

where w q = v^lql, |q| <C 1, is the magnon energy, and 
P n , P a are the normal and anomalous polarization opera- 
tors which can be calculated perturbatively for the inter- 
action Eq. JSJ. The above Green's function is the normal 
one (the anomalous Green's function D a (u>,q) has the 
same denominator). 

The magnetic stability requires that all the poles of 
(|19|l are at positive uj 2 , i.e. no imaginary poles exist. 
The one loop, first order contributions to P n ,P a were 
already calculated in Ref. fl3L and we now calculate the 
second order, to match the order in the ground state en- 
ergy. Due to the fact that in the spiral state the fermion 
operators on the two sublattices are mixed via Eq. (7J, 
corrections to the vertex appear in lowest order. For 
example the second order diagrams for the normal part 
Pn are shown in Fig. 5(a,b,c). After evaluating these 
diagrams we obtain: 



Mi 



i^)(0,q)^ £ 



+ - 




(™k - «k+q)(«p - «p+q) 



p,k,AT p \ ^k e k+qA e P £ p+qi 



(^kVqX^+^p) 

2n k n p 

(4-^X4-^), 



(20) 



For the stability analysis one needs only the long wave- 
length, i.e. |q| — * limit of (|20|l . In this limit the second 
and third terms (Fig. 5(b,c)) in (|2(Jfl cancel out, leaving 
only the first term, Fig. 5(a). The corresponding diagram 
for the anomalous part (Fig. 5(d)) is: 



Pi 2 )(0,q) = -e- 2 ^i 2 >(0,q), |q| 



0. 



(21) 



Adding to these results the one- loop results, P^l from 
Ref. IT^L we obtain the polarization operators up to two 



loops P nt a — Pn)a + Pn,L- For |q| — > we have explicitly: 



,(2) 



Pn(0,q) = 



V2Z 2 t 2 



2e F 



7T 2 /?l/3 2 



|q|, 

(22) 



(a) 



(b) 



(c) 





(d) 




+ 



FIG. 5: (a,b,c) Two-loop, second order diagrams contributing 
to the normal polarization operator P„. (d) Two- loop contri- 
butions to the anomalous polarization P a (only first diagram 
of several is explicitly shown). 



D /n \2i» V2Z 2 t 2 ( 2e F \. 4V2IZH 4 
P a (0,q)e 2 ^ = — ? == ( 1- — ] |q| - - |q| 



The quantity I is defined as: 



7T 2 0lP2 



(23) 



' dxdy 



(cosx — cosy) 2 



o (27r) 2 (cosx — cosy) 2 + (/3i//32)(sinx— s'my) 2 
The stability criterion reads: 

[^ 2 +2o; q P„(0,q)] 2 >4 W 2 |P Q (0,q)| 2 , (25) 

which, for the (1,0) spiral state, is equivalent to the fol- 
lowing inequality, in terms of the microscopic parameters: 



2Z 2 t 2 
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7r 2 /?l/? 2 

8IZH A 



TT 2 PlP2 



(26) 



We have verified numerically that inside the region 
marked "Stable" in Fig. 4, the expression inside the ab- 
solute value sign on the left hand side of (|2l))l is positive, 
whereas the one on the right is negative. This means 
that (|26|l is equivalent to the equation 



1 - 4p, > . 



(27) 
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This is the same condition as the one obtained^ in the 
one-loop approximation. This condition is always sat- 
isfied since 1 — 4p s = 1 — Z p , and Z p = 0.72 due to 
spin quantum fluctuations renormalization in the un- 
doped S=l/2 antiferromagnet. The physical meaning 
of the "quantum stabilization" is quite simple. Indeed, 
treating the spins semiclassically we would have p a = 1/4 
and hence a marginal ground stateAi2ii£ The spin quan- 
tum fluctuations reduce the spin stiffness and hence make 
the energy of the spiral state lower, see Eq. 0. This 
is certainly only the intuitive physical picture, while 
the regular proof is the one presented above, leading to 
Eqs. ESEJ . 



IV. SELF-ENERGY CORRECTIONS AND 
THEIR INFLUENCE ON STABILITY 

The last issue we would like to address is the variation 
of the hole dispersion with doping and its influence on the 
phase diagram of Fig. 4. While the renormalization of the 
one-hole properties at zero doping has already been taken 
into account via the SCBA, Eq. (J3J, at finite doping the 
many-body hole-hole interaction corrections originating 
from the vertex H10|) have to be calculated. 

The leading self-energy diagram, contributing to the 
hole dispersion at finite doping is shown in Fig. 6. The 
corresponding expression is 



8Z 2 t 2 



2-^ it 



(28) 



The self-energy is the same for tp and <p fermions and 
therefore it influences the dispersion 0, but does not 
influence the gap A between the bands. The one-loop 
self-energy depends only on the momentum of the parti- 
cle. We remind the reader that the one-loop self-energy 
is closely related to the one-loop correction to the total 
energy, namely by cutting a fermionic line in diagram 
Fig- da) one obtains the self-energy Fig. 6. This means 
that the self-energy effect has already been taken into ac- 
count in the calculation of the total energy performed in 
Section II. Nevertheless it is interesting to see explicitly 
how the dispersion changes. The integral in 1|28|) can be 
easily evaluated numerically at any point of the phase di- 
agram Fig. 0] However it is more instructive to consider 
the case of isotropic dispersion /3± = (3% = j3 where the 
integration can be performed analytically. In this case 



n,. 



2Z 2 t 2 e F 
tt/3 

f 1 2 



[i + n k {kl-k 2 2 )} 



i f-i 

W \ x 2k 7 J 



for k < 
for k > 



kp 
kp 



(29) 



With account of the self-energy the dispersion should 
be replaced to — > + E^. Thus effectively we have a 
deformation of the Fermi surface 

7 2 t 2 7 2 t 2 
f3i=P+ , & = P • (30) 




FIG. 6: One-loop self-energy correction. Dashed line repre- 
sents the vertex . The correction is the same for -0 and if 
fermions. 



Notice that the deformation does not depend on the 
doping S. On the other hand numerically this effect turns 
out to be rather weak. For the purpose of an estimate we 
take t' = —0.8, t" = 0.7, corresponding to the physical 
values of these parametersi 2 ^ 2 ^ At this point fix = 2.95, 
fa = 3.80, Zt = 1.19, hence the Fermi surface deforma- 
tion is Ap//3 w Z 2 t 2 /{ny/pTP2) « 0.13. 

Finally, we can summarize our findings: (i.) The lead- 
ing order self-energy correction produces a weak, doping 
independent spectrum renormalization. (ii.) The charge 
stability analysis of Section II performed up to (two-loop) 
order (Zt) 4 is not influenced, at that order, by the change 
in the spectrum. Taking the latter effect into account 
produces higher order diagrams and is expected to be a 
very small correction, due to the previous point (i). (iii.) 
The magnon stability analysis of Section III is not influ- 
enced by the inclusion of l|29|l . since the the density of 
states cx l/^/PTJh cancels out in the stability criterion 

G3- 



V. CONCLUSIONS 

We have studied the stability of the spiral states in the 
t — t! — t" — J model at low doping S <C 1. The stability 
was studied both in the charge and in the spin sectors, 
and takes into account interaction effects by including 
diagrams up to two loops. Our main result is that rela- 
tively small values of t',t" are sufficient to stabilize the 
(1,0) spiral state, as shown in Fig. 4. For example for 
t" = 0, the critical value of t' is \t' c /J\ ~ 0.25, or in 
units of t, \t' c \ w 0.08i (for t/J = 3.1). The unstable 
region for small t' , t" is a good candidate for an inhomo- 
geneous ground state of some sort, e.g. a striped one. 
Our results are also consistent with DMRG studies^ of 
the extended t — J model predicting a transition from 
stripes at t' = to a homogeneous ground state around 
\t' /t\ w 0.1. Parameters for real cuprates correspond to 
the top left corner of the phase diagram Fig. 4 where the 
(1,0) spiral state is stable. 

We have shown that the interaction effects (higher loop 
diagrams) should be taken into account and they tend 
to increase the region of spiral stability in the space of 
parameters t',t". The effective coupling constant gov- 
erning the perturbation theory is g — (Zt) 2 /it sa 0.3. 
Therefore each successive order is several times weaker 
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than the previous one, and thus we expect our results, 
calculated up to second order, to be quite reliable. The 
good convergence is demonstrated by the inset in Fig. 4 
where we show the critical lines calculated in zero, one 
and two-loop orders. The magnetic stability of the spiral 
state is guaranteed by spin quantum fluctuations (order 
from disorder effect) - this conclusion remains valid in 
the two-loop approximation. 

Our calculations are based on the chiral perturbation 
theory, therefore they are reliable only in the limit 6 -C f , 
whereas the finite doping renormalizations of the various 



Green's functions and vertices become more and more 
substantial as doping increases, leading essentially to an 
untreatable problem. 
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